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INTRODUCTION 
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The release of volatile organic compounds (VOC's)into 
the atmosphere has come under fire by the Environmental 
Protection Agency (EPA) and state regulatory agencies in the 
past several years. These releases come from a wide variety 
of sources from gasoline filling stations to industrial and 
military operations. There are also countless natural 
sources of VOC's including volcanos and forest fires. 
Environmental groups have had a great influence on 
government agencies resulting in the enactment of numerous 
statutes and regulations to reduce these releases over an 
approximate five (5) year period. 

Enforcement of environmental regulations is carried out 
by the individual states which impose their own specific 
regulations which in some cases are even more stringent than 
those of the federal government. A case in point is the 
state of California where regulations are resulting in rapid 
actions on the part of corporations to comply or face harsh 
penalties. These actions include both a switch to 
alternative manufacturing materials and employment of new 
pollution control devices. 

A new pollution control technology currently under 
development utilizes ultraviolet (UV) light to break down 
the volatile organic compounds (VOC's) contained in an 
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airstream. One pollution control equipment manufacturer has 
already employed this method in an installation in 
California (see Figure (1) ) . The process involves injecting 
ozone into a mixture of gaseous VOC's which have been 
radiated by UV light. In the presence of oxidants and UV 
light, the air/VOC mixture undergoes a series of photo- 
dissociation and radical oxidation reactions which destroy 
the VOC's. Despite this pollution control equipment's 
success in controlling a wide range of airborne VOC's, the 
extent of VOC reduction due to photolysis in the photolytic 
reactor section of the system is unknown. Our lack of 
information on the effectiveness of this photolytic reactor 
as shown in Figure (2) is due to proprietary information 
retained by the manufacturer, insufficient test data, and 
lack of an adequate model to predict its performance. 

There have been some analytical and experimental 
studies reported in the literature on photochemical reactors 
(ie., Y. Harano and T. Matsura (1972)). However, these 
studies are generally directed towards simplified geometries 
which can be modeled analytically (ie., axial, elliptical, 
and columnar reactors) and which assume well-mixed 
conditions, low light absorption, uniform flow, and 
negligible diffusion. Experimentation has been directed 
towards determining the deviation between these highly 
simplified models and actual reactor performance in order to 
determine the effects of factors omitted. Few if any 
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TERR-AQUA ENVIRO SYSTEMS 
AIR POLLUTION CONTROL TECHNOLOGY 



TAES ULTRAVIOLET-OXIDATION CONTROL SYSTEMS have been designated as Best 
Available Control Technology (BACT) for Volatile Organic Compounds (VOC) and Reactive Organic 
Compounds (ROC). 

The UV-Oxidation processes utilize a proprietary, on-site, self-regenerative system designed to meet 
the current EPA and most stringent local Air Quality District rules for capture and destruct 
efficiencies. 



TAES AIR POLLUTION CONTROL FLOW DIAGRAM 




Figure (1) - VOC Pollution Control Device Schematic 

(Terra-Aqua Environmental Systems (1992) ) 










Figure (2) - VOC Pollution Control Device Photolytic Reactor 



5 



studies can be found in the literature which attempt to 
quantitatively model photochemical reactors of complex 
geometries and flow fields using numerical methods. This 
report provides a numerical model which can be used to 
quantitatively describe the photo-dissociation/radical 
oxidation of VOC' s in the photolytic reactor of Figure (2). 

A secondary but equally important objective for this 
analysis is to develop an analytical tool which can be used 
to optimize photo-reactor geometry. This numerical model, 
once verified with experimental data (in progress at Penn 
State) and expanded to incorporate additional reactants, can 
be utilized to predict optimum geometry for the photolytic 
reactor section of a pollution control device. 

Finally, this model will assist the designer in 
determining when the use of a photolytic reactor in a 
pollution control device is appropriate. Since there are a 
multitude of VOC' s which vary widely in their ability to 
react when exposed to UV light, cost effective 
implementation of a photolytic reactor may be highly 
dependant on the specific VOC's expected in an airstream. 

For example, the photolytic reactor may not provide any 
additional benefit in a pollution control device which 
treats specific emissions from a military or industrial 
painting operation. This numerical model may indicate when 
a significant cost savings is possible by omitting the 
photolytic reactor section and could also predict when the 
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photolytic reactor can function as the sole VOC reduction 
device in a pollution control system. 

1 . 1 Literature Review 

A review of literature on the subject of design and 
analysis of photochemical reactors shows that very little 
work has been accomplished in the area of numerical 
modeling. 

Jacob, S. M. and Dranoff, J. S. (1969) conducted 
experiments on light distribution within an elliptical 
photoreactor and compared these results with some closed 
form solutions. They found that the closed form solutions 
are unable to accurately predict light intensity within an 
elliptical reactor. Additional studies by Jacob, S. M. and 
Dranoff, J. S. (1970) concluded that numerical means are 
required to accurately determine the light intensity 
distribution within a well-mixed elliptical photoreactor. 
While pointing out the necessity to use numerical methods in 
the analysis of photolytic reactors, neither of the above 
studies attempted to deal with reactors containing non- 
uniform flow fields and spatially varying reactant 
concentrations. In both studies only the case of simple 
elliptical geometry with laminar flow and well-mixed 
conditions was considered. 

Work by Ragonese, F. P. and Williams, J. A. (1971) on 
the scale-up of a laboratory photoreactor model to meet 
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actual chemical production requirements outlines the 
difficulties in utilizing both dimensional analysis and 
empirical data. Although success has been achieved in this 
scale-up process, the methods used apply to simplified 
geometries with laminar flow and well-mixed conditions. 

This study points out the need for numerical modeling to 
enable a designer to scale photoreactors to meet production 
needs or to achieve sufficient pollutant reductions in a 
pollution control device. 

The literature review provided by Harano, Y., and 
Matsura, T. (1972) describes the many difficulties involved 
in the design of photochemical reactors. They summarized 
that even for the simplified geometries which have been 
studied extensively to date including cylindrical, parallel 
plate, annular and elliptical reactors, the actual 
concentration gradients, light intensity distributions and 
flow conditions greatly influence overall reaction rates 
within the photoreactor. 

All the literature reviewed on the subject of 
photolytic reactor design points to the need for numerical 
reactor modeling which incorporates both photochemical and 
fluid mechanics principles and the latest numerical solution 
techniques . 



Chapter 2 

PRELIMINARY MODELING ASSUMPTIONS 
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Due to the great difficulty in modeling the entire 
reactor under investigation including all possible chemical 
reactions, this model describes only the photo- 
dissociation/radical oxidation (photolysis) of formaldehyde 
in a rectangular flow field perpendicular to the axis of a 
single UV lamp as shown in Figure (2) . Formaldehyde has 
been chosen for this model since its photochemistry has been 
studied, documented and verified. Although this model 
applies only to the photolysis of formaldehyde, it is 
general enough to incorporate the chemistry of other 
airborne components as well as combinations of other 
photochemically reactive solvents in the future. In 
addition, the principles required to expand this numerical 
analysis to additional materials and/or chemical reactions 
have been included in this report. 

The following assumptions have been made to simplify 
the analysis: 

A. Constant fluid properties (ie ., density , 
temperature, diffusion coefficient) . 

B. Negligible thermal effects (ie ., solution of energy 
conservation equations not required) . 

C. Two-dimensional flow field. 
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D. Weak absorbance of LJV light (ie., light intensity 
is independent of formaldehyde concentration) . 

E. Photo-dissociation/radical oxidation of 
formaldehyde only. No other reactions between formaldehyde 
and other compounds are included except those required to 
complete the photo-dissociation/radical oxidation of 
formaldehyde . 

F. Steady, time-averaged flow field (ie, no 
oscillation of vortex shedding in the wake region) . 

G. No wall losses (ie., chemical reactions which 
sometimes occur due to collisions of molecules into the 
reactor walls or the UV bulb itself) . 

H. No reflection of UV light from the reactor walls. 

The above assumptions are considered to be reasonable 
ones for use in this initial numerical model. In the 
future, factors such as wall losses and reflection can be 
introduced into the numerical simulation to increase its 



accuracy. 



Chapter 3 

MASS CONSERVATION WITHIN THE PHOTOLYTIC REACTOR 



10 



Given the assumptions provided in Chapter 2, the 
differential equation which describes material concentration 
in the photolytic reactor can be written: 

-v°Vc + DV 2 C + n = dC/dt [1] 

where, C = concentration of formaldehyde (PPM) . 

v = flow velocity (with components v x ,v y ) (m/s) . 

D = formaldehyde diffusion coefficient (m 2 /s). 

Q = generation of formaldehyde (PPM/s). 

The first element of Equation [1], the convection term, 
involves the fluid mechanics of the reactor. It combines 
concentration changes along a given direction with the fluid 
velocity to result in a time rate of change in 
concentration . 

The second term in the above equation is the diffusion 
term. It describes the time rate of change of concentration 
entering/leaving the control volume due to concentration 
gradients at the element surface. Note that this term is 
dependant on the second spatial derivative of C multiplied 
by a diffusion constant and is relatively small with respect 
to the velocities in the flow. To achieve the highest 
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accuracy possible for all possible concentration gradients 
and flow velocities, this term is included in the model. 

The third term in this equation describes concentration 
changes due to the photolytic action of the UV bulb and 
other chemical reactions. For this model the "generation 
term" simplifies to the following: 

Q = -KC [2] 

where, K is the overall specific absorption rate constant 
(overall rate constant for the breakdown of formaldehyde by 
UV light) . 

The final term is the resultant time rate of change for 
the element which, when included yields a transient 
solution. Since the transient solution has practical value 
both in predicting the start-up concentrations within the 
reactor as well as in examining time varying inputs, the 
time dependant term has been included. In addition to the 
above benefits, a transient solution must be found in order 
to obtain convergence for the numerical solution technique 
used. This is required since it is difficult for the 
iterative numerical technique (Gauss-Seidel) to jump 
directly from zero time to an infinite solution (steady- 
state solution) without some intermediate steps in time. 

For the case of a multi-component flow field which 
includes chemical reactions as well as photolytic reactions. 
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Equations [1] and [2] can be written: 

-voVC ± + D-jV 2 ^ + Q. = dCi/dt [3] 

+ other reactions (see Appendix B) [4] 

Analysis of the multi-component case requires the solution 
of n simultaneous Equations ([3] and [4]) corresponding to 
the total number of reactants. Note that Appendix B 
contains numerous additional equations involving chemical 
and photochemical reactions between many possible airborne 
constituents. In order to simplify this analysis only the 
photolysis of formaldehyde will be considered (ie., no 
reactions between formaldehyde and other airborne 
constituents will be considered) . 

In the next chapter the determination of velocities 
required in the convection term of Equation [1] is 



discussed. 



Chapter 4 

FLUID MECHANICS PRINCIPLES 
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For the case of turbulent, viscous flow around the 
spherical bulbs shown in Figure (2) , a complete solution of 
the Navier-Stokes equations is required to obtain accurate 
fluid velocities. Due to the difficulty in solving the 
complete Navier-Stokes equations, a commercially available 
computational fluid dynamics program has been employed to 
determine the velocities in a rectangular region around a 
single bulb. Figure (3) shows a standard K-£ method solution 
of the Navier-Stokes equations for flow of air around a 
circular cylinder generated by the FLUENT program. The 
shape of a 0.03m diameter circular cylinder in a 0.05m x 
0.20m flow field has been approximated by a series of grid 
points (100 wide by 25 high) . The program calculates x and 
y velocities at every grid point in the flow for input into 
the convection terms of Equation [1] . 

The following parameters were input to the FLUENT 
program (properties of air/solution parameters) : 

Properties of air only assumed (no air/form, mixture) 

Length of domain (x direction, m) - 0.2 

Height of domain (y direction, m)- 0.05 

Diameter of Bulb (m) - 0.03 

Grid size (x dir, y dir)- 100 x 25 



14 




Figure (3) - FLUENT Generated Streamlines For Flow of 



Air Around a UV Lamp in a Uniform Flow Field 
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Reynolds Number, Re 0 - 4,500 

Inlet Temperature (°C)- 35 

Inlet Pressure (Kpa)- 101.33 

Inlet Density (kg/m 3 )- 1.293 

Inlet Viscosity (kg/m-s)- 1.72 x 10 -5 

Inlet Turbulence Intensity (%)- 10.0 

Inlet Velocity (m/s)- 2.0 

Since it is beyond the scope of this report to examine 
the detailed accuracy of the FLUENT program output, these 
results will be assumed to be sufficiently accurate for the 
purpose of making design decisions. There are a variety of 
models and techniques which can be employed to increase and 
verify the accuracy of the flow analysis. However, since it 
is the purpose of this report to model the photolytic 
reactor in a fashion which well-exceeds the accuracy of 
simple plug flow models which can be found in the 
literature, the unmodified FLUENT output is used. 



Chapter 5 

CHEMICAL PRINCIPLES 
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The chemistry associated with the photolytic reactor 
under consideration here is much the same as that of 
atmospheric chemistry with the exception that energy is 
added with a mercury lamp which emits several discrete 
wavelengths (predominantly 253nm) versus sunlight which 
contains these and other wavelengths. The chemistry of the 
photolytic reactor can be broken down into two categories 
including photochemistry and free radical chemistry. 

5 . 1 PHOTOCHEMISTRY 

Photochemical processes initiated by UV radiation begin 
with the absorption of a photon by a molecule as follows: 

K 

A + hD — » A* [ 5 ] 

where, A* is an exited A molecule and hu is a photon of UV 
energy. Upon excitation the A* molecule can undergo the 
following reactions: 



K 1 

Dissociation: A* — > B + C [6] 

(A* broken down into stable molecules by itself) 

K 2 

Direct Reaction: A* + B — » D + E [7] 



(A* combines with B to form two stable molecules) 
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K 3 

Fluorescence: A* — >A + hi) [8] 

(A* returns to orig. state after release of energy) 

Collisional K 4 

Deactivation: A* + B — >A + B + hi) [9] 

(A* returns to orig. state after collision with B) 

Each one of the above reactions has a probability of 
occurrence called a quantum yield, . For example, a 
quantum yield of ^ = 0.5 implies that there is 50% chance 
that A* molecules will undergo dissociation (Equation [6]) 
while there is a 50% chance A* molecules will directly 
react, undergo fluorescence or deactivate from a collision. 
Since all the A* molecules must undergo one of the above 
reactions : 

h 

E <!>• = 1 [10] 

hi 1 

The production of A* molecules is given by the 
following rate equation: 

d fA*1 = K [A] [11] 

dt 



In Equation [5], since one A* molecule consumes one 
molecule of A: 



d [A] = -K [A] 
dt 



[ 12 ] 
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The symbol K is the overall "specific absorption rate" 
constant (also given in Equation [2]) and [A*], [A] are the 

concentrations of A* and A, respectively. Note that the 
overall rate constant corresponds to all possible quantum 
yields for A*. By combining equations [5] through [9], the 
rate equations corresponding to the dissociation, direct 
reaction, fluorescence, and collisional deactivation 
reactions (termed the first order rate equations) can be 
written as follows: 



d r A 1 = -4> 1 K [A] = -K-jJA] (reduction in A due to [13] 

dt dissociation) 

d fA] = -<J> 2 K[A] = -K 2 [A] (reduction in A due to [14] 

dt direct reaction) 

d f A] = -<E> 3 K[A] = -K 3 [A] (reduction in A due to [15] 

dt fluorescence) 

d [A] = -<I> 4 K[A] = -K 4 [A] (reduction in A due to [16] 

dt collisional deactivation) 



The quantities K 1 (=0-^K) , K 2f K 3 , and K 4 in Equations [13] 
through [16] are termed the first order rate constants. 

Note that multiplying both sides of Equation [10] by K 
results in K]_ + K 2 + K 3 + K 4 = K. Therefore, the overall 
reaction rate for the reduction of A is equal to the sum of 
the individual rates for reactions which can occur upon 
excitation of A molecules. 

For the general case of atmospheric chemistry, the 



first order rate constants can be determined through the 
following equation: 
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K 1 = \ a(A.,T)<D^ (A.,T) I (A.) dA. 



[ 17 ] 



a(A.,T) = absorption cross section of the A molecule at 
operating temp., pressure and UV wavelength 
(cm 2 ) . 

^(XjT) = quantum yield or probability that molecule A will 
dissociate upon excitation (recalling that one 
molecule of A produces one molecule of A*) (K2 
would correspond to direct reaction rate) . 



I (A.) dA. = irradiance of incident UV light (photons/cm 2 -sec) . 



The above first order rate constant, K-^ corresponds to 
the rate at which molecule A will undergo dissociation (one 
of the four reactions mentioned above) when exposed to a 
broad spectrum of light energy in the atmosphere. 

In order to calculate the first order rate constants 
for combinations of specific, discrete wavelengths the 
following approximation is utilized: 

n 

K x = E a(A. i ,T)0) 1 (A^T) I (A. i ) AA. 



[ 18 ] 
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The quantities a and ^ above are average values at . The 
quantity I(X i )AX is sometimes expressed as J(A.^) (actinic 
flux) which results in the following: 

n 

k x = E ca^T^a^T) j(X ± ) [i9] 

C = l 

5.1.1 DETERMINATION OF FORMALDEHYDE PHOTOLYTIC RATE 
CONSTANTS 



Formaldehyde participates in the following photolytic 
reactions (A.<300nm) : 



K 1 

HCHO + hv — »H + HCO (HCHO is formaldehyde) [20] 



K- 



H 



2 + CO 



[ 21 ] 



The wavelength of UV light utilized in this device is 
predominately X = 253 nm. The bulb dimensions are as 
follows : 



D (diameter) = 0.033 m 
L (length) = 1.5536 m 

Ptotal ( tota i bulb power requirement) = 40 W (J/s) 
^visible (bulk output in visible spectrum per 
manufacturer) = 40% Ptotal 

p 253 (bulb output at 253 nm) = 95% P v i s ibi e 
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As an example, the first order rate constant at ^.=253nm 
(corresponding to the reaction of equation [20] above) is 
determined as follows: 

E = hc/A. = 6 . 625 6x1 0~ 34 J-s /photon (2 . 997 9xl0 8 m/s) 

253x10 

= 7.8509 x lO'^J/photon 

A BU lb = 7tDL = Jt(33mm) (1553.6mm) = 1611 cm 2 
P 253 = 40W (40%) (95%)= 14.44 J/s 

J(A=253nm) = _P 2 53 = 14.44 J/s 

Abulb (E) 1611 cm 2 (7.8509 x 10~^J/photon) 

= 1.1419 x 10 16 photons/cm 2 s 

a(253nm,25C) = 3.79 x 10~ 21 cm 2 /molecule (R. Atkinson, 

et al . (1992) ) 

<$! (253nm, 25C) = 0.3 (R. Atkinson, et al . (1992)) 

= G(X,T)0 1 (X,T) J(A) 

= 3 . 79xl0 _21 cm 2 /molecule (0 . 3) (1 . 1419xl0 16 phot . /cm 2 s) 

= 1.30 x 10~ 5 (1/s) 

K = overall specific absorption rate (yields both 
reactions 1,2) 
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K = K x + K 2 = a (A., T)0 -l (A,, T) J (A,) + G(l, T)<I>2 (X, T) J(X) 

= a(X,T) (O l a,T)+d> 2 (^,T) ) J(X) 

<I >2 (253nm, 25C) = 0.5 (R. Atkinson, et al . (1992)) 

K = (3 . 79xl0 -21 cm 2 /molecule) (0.3+0. 5) x 

(1 . 141 9xl0 16 photons/cm 2 s) 

= 3.47 x 10 -5 (1/s) (yields both reactions 1,2 

(Eqns [20] , [21] ) ) 

5.1.2 LIGHT INTENSITY DISTRIBUTION THROUGHOUT THE REACTOR 

The rate calculations above are based on the intensity 
at the bulb surface. This formulation is based on the Beer- 
Lambert law applied to weak absorbance conditions such as 
the atmosphere (Finlayson-Pitts, B. J. , and Pitts, J. N., 

Jr. (1986) ) . Since the energy added via the bulb per unit 
area and corresponding rate constants vary indirectly with 
the radial distance from the bulb surface. Equation [2] must 
be modified as follows: 

Q. = -K(R/r)C [22] 

where, r = radial distance from bulb centerline. 

R = radius of bulb. 

In addition, the light intensity in an actual 
photochemical reactor is also attenuated according to the 
well known Beer-Lambert law (Finlayson-Pitts, B. J. , and 
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Pitts, J. N., Jr. (1986)) which is written as follows: 

I = Io exp{-aNP } [23] 

where, Io = monochromatic source intensity (photons/cm 2 -sec) . 
o = absorption cross section (cm 2 /molecule) . 

N = molecules/cm 3 (molar concentration) . 

P = path length from the source (cm) . 

For low concentrations and/or low absorbance, as is the 
case with this model, the effects of concentration on light 
intensity can be shown to be negligible. For this reactor: 

a = 3.79 x 10 -21 cm 2 /molecule 

N = 1.229 x 10 15 molecules/cm 3 (equivalent to 50 PPM) 

P = approx. 20 cm 



Substituting into Equation [23] yields: 



I 



Io expl-(3.79 x 10 21 cm 2 /molecule) x 

(1.229 x 10 15 molecules/cm 3 ) (20 cm)} 
0.99991 Io 



In view of the above result it is clear that the weak 
absorbance assumption is a reasonable one. 
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In order to accurately model a multi-component air/VOC 
mixture where light absorbance must be considered, the 
following multi-reactant Beer-Lambert law is utilized: 

I = Io exp{-[a 1 N 1 + C 2 N 2 + ] L } [24] 

Using this relationship an intensity function, I = f(x,y) 
can be determined using the actual concentration of 
materials in the reactor which are capable of absorbing 
light . 

The specific photochemistry of reactants which could be 
considered in this model are contained in Appendix B. 
However, for the purpose of this initial model only the 
photolysis of formaldehyde is under consideration. 

Therefore, the intensity will be considered to be a function 
of the radial distance from the bulb surface only due to the 
weak absorbance of formaldehyde which was demonstrated 
above . 

5.2 FREE RADICAL CHEMISTRY 

As can be seen in Appendix B there are a multitude of 
reactions which can be included in this model as there are 
numerous materials normally contained in an airstream which 
are capable of undergoing photolysis and reacting with 
formaldehyde. However, due to the great complexity of 
including all these reactions, this initial model will be 
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only concerned with the removal of formaldehyde through 
photolysis alone. Therefore the "generation" term of 
Equation [1] will be that given by Equation [20] with no 
other reactions included. 

Note that in making the above simplification, reactions 
with free radicals have been ignored. Therefore, the 
results of this investigation apply only to photolytic 
reactions. This analysis is a logical first step toward a 
more comprehensive analysis that will include free radical 
reactions. This analysis also provides valuable insight 
into the photolysis of formaldehyde in a non-uniform flow 
field. 



Chapter 6 

CLOSED FORM SOLUTION FOR SIMPLIFIED FLOW FIELD 
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In order to get a general idea of the size of the 
actual pollution control device required to reduce a given 
formaldehyde input concentration by a factor of two through 
photolysis, a simplified flow field model can be employed. 
Restating the mass conservation equation: 

-v°Vc + DV 2 C + Q = dC/dt [ 25 ] 



The following simplifying assumptions will be made: 

A. Two-dimensional, inviscid, steady flow. 

(v = velocity in x direction = constant, 
velocity in y direction = 0) 

B. Constant radiance intensity throughout the reactor. 

C. Formaldehyde photolysis only. 

D. Diffusion negligible. 

With these assumptions, the above equation simplifies to: 

-v°VC + Q = 0 [26] 

Since the only reaction involves the photolysis of 
formaldehyde, substitution of Equation [2] into Equation 
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[26], performing the dot product and simplifying yields: 



v dC = -KC [27] 

dx 



Rearranging terms and solving for C(x) yields: 



In (C) -In (C Q ) = -Kx/v 


[28] 


ln(C/C Q ) = -Kx/v 


[29] 


C/C Q = exp{-Kx/v] 


[30] 



The values assumed are: 

C/C Q =0.5 

K = 3.47 x 10 5 (1/s) (as calculated) 
v = 5 ft/s 

Substituting in values and solving for x yields: 

Xcjq = 99, 877 ft = x corresponding to a 50% C reduction 
= 18.9 miles 

This result is a very pessimistic one which predicts 
poor performance for a reactor which relies on the 
photolysis of formaldehyde alone. However, in examining the 
general solution of the analytical model for different 
values of V and K, it appears that the analytical model can 



Fornt*1d*hy de Concentration (PPM) 
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Figure (4) - Formaldehyde Concentrations in an Idealized 



Two-dimensional Columnar Photo-reactor 



be significantly improved (see Figure (4)) by employing a 
higher intensity mercury lamp (increasing K) and/or by 
decreasing the flow velocity. 

Increasing the intensity of the lamp is a viable 
option. There are currently bulbs available which are 
capable of over 1000 times the power of the mercury vapor 
lamps which are used in the photolytic reactor of Figure 
( 2 ) . 

Since a decrease in velocity beyond some practical 
value will only lead to a decrease in the amount air 
treated, increasing the bulb intensity is clearly the best 
practical alternative. 
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Chapter 7 

DEVELOPMENT OF THE NUMERICAL MODEL 



The numerical photolytic reactor model under 
consideration predicts the concentrations in a rectangular 
area around a typical bulb in the reactor of Figure (2) . 
Since very little reduction in formaldehyde is expected in a 
small 21.2 x 6.6 cm domain utilizing the actual rate 
constant calculated in chapter 5.1.1, the rate constant has 
been increased 10,000 times to an attainable 0.35 (1/s) to 
achieve some discernable results within two decimal places. 
Although this bulb intensity/reaction rate increase is 
presently impractical from a cost standpoint (very expensive 
bulbs are required) it greatly improves the analysis and 
interpretation of the numerical results for the small flow 
domain under consideration. The following data describe the 
parameters of the model: 

L, length of the domain in x direction = .2121 m 

M, number of nodes in the x direction = 100 

h, size of increments in the x direction = 0.002143 m 

H, length of the domain in the y direction = .0664 m 

N, number of nodes in the y direction = 31 

h, size of increments in the y direction = 0.002143 m 

Uo, inlet velocity = 2.0 m/s (const, across entire 



area) 
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Co, inlet formaldehyde concentration = 50.0 PPM ( " ) 
Bulb diameter = .03 m 

K, overall rate constant at bulb = .35 1/s 
D, diffusion coefficient = 0.00002 m 2 /s 

£, solution convergence criteria = 0.000001 PPM 

(maximum difference between successive concentration 
calculations at each point at a given time) 

At, time increment = 0.01 s 



The computer model uses a finite difference technique 
in conjunction with a Gauss-Seidel iterative method to 
calculate the concentration field throughout the reactor at 
selected time intervals. The mass conservation equation is 
restated as follows: 

-voVc + DV 2 C + Q = dC/dt [31] 



Since the convection terms dominate Equation [31], an 
upwind difference method has been incorporated into the 
above equation to improve the accuracy of the results. 

Simply stated, the upwind technique allows the concentration 
at any point i,j to be calculated based on points upwind 
since concentration changes are much more closely related to 
the fluid flow than they are to diffusion and photolytic 
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effects. Hence, differencing in the upwind direction is 
utilized for the convection terms. For the sake of 
simplicity, central differencing is used for the diffusion 
terms. The convection term of Equation [31] will now be 
developed in terms of finite differences. 

The equation which follows, Equation [32], is an upwind 
difference expression which is valid when both the fluid 
velocity components, v x and v y are positive. 

-v°Vc - -<v x /h> !> [32] 

When both the fluid velocity components are negative: 

-V°VC = - (v x /h) (C i+lf j-C if j)-(v y /h) (C ±/ j+1 -C ±/ j) [33] 

An expression derived from Equations [32] and [33] above 
which will always yield an upwind difference for any 
combination of velocity component directions can be written 
as follows: 



-voVC = 




- (l/2h) I <|v x |+ v x ) 

|+ v y ) (c if j -C i/ ;)_!> + (V x - |v x |)(C i+1/ 
(v y - | V y I ) (C i# j+1 -C i/ j)l 





+ 



[34] 



While the above equation may appear unnecessarily lengthy. 
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it greatly increases the speed of numerical calculations as 
a decision does not have to be made (FORTRAN IF statement) 
at each point in order to apply the correct equations based 
on the velocity component directions. Without condition 
statements in the calculation of the convective terms, much 
of the program can be vectorized to greatly reduce the 
computation time. 

In terms of central differences, the diffusion term can 
be written: 



DV C (D/h 2 ) I C i+1 ^ j + C^-^j 



^i, j + ^i,j+l + *^i, j-lJ 



The remaining terms are: 
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Q. = -KC.j , [36] 

± e J 

dC/dt = (1/At) (C ± , n - Ci i n_1 ) [37] 

1 r J 1 1 J 

The superscripts n and n-1 above refer to the present 
value and the value at the previous time, respectively. It 
is important to note that a backward difference scheme has 
been employed for the time dependant solution of the 
governing equation. Although all terms on the left side of 
the governing equation (Eqn. [31]) have been developed using 
concentration values at the present time, the superscript n 
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has been omitted for clarity. The backward difference 
method results in an implicit set of equations which can be 
solved by the Gauss-Seidel iterative method to find the 
concentration at each point in the reactor at a particular 
point (iteration) in time, n. This method may seem 
cumbersome in comparison to the forward difference method 
which can explicitly calculate the concentrations at a 
future point in time without iteration. However, the 
backward difference method has higher stability which allows 
a much larger time step to be used. 

Substituting Equations [33] through [37] into Equation 
[28] and solving for j yields the following iterative 
equation : 



The superscripts k and k+1 in equation [38] represent 
concentration values calculated at the past and present 
iterations, respectively. Note that the C i _ 1 ^j and 
terms have the superscript k+1 whereas the C^ +1 ^ j and 
Ci,j + i terms have the superscript k. This is due the 
iterative "search pattern" which sweeps from j=l to 




1 



(1/2)1 ( I v x | + v x ) 



|v x |+ |v y |+ (4D/h) + hK + (h/At) 




[38] 
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N (y=0 to y=H) for each i until i=M (x=L) is reached. This 
particular sweep pattern in conjunction with the Gauss- 
Seidel method makes use of the latest values (k+1) for 
C^-^j and in calculating the updated C i ^j k+1 value. 

The boundary conditions of the flow domain are set as 
follows : 



A. Inlet - C i = const. = Co (i=l;l<j>25) [39] 

B. Top and bottom surfaces - dC/dy = 0 (symmetry) or, 

c i,j+l = c i,j-l (1— i-100; j=l, 25) [40] 

C. At bulb - dC/dx, dC/dy = 0 (solid boundary) or, 

C i+1 ,j = (see program, Appendix A) [41] 

(see program. Appendix A) [42] 

D. Exit - dC/dx = 0 (minimal change at exit) or, 

c i,j = c i-l,j (i=100;l<j>25) [43] 



The computer program included in Appendix A uses 
Equation [38] (including the relationships of Equations [39] 
through [43] when required) to calculate and update new 
concentration values at all grid points in the domain per 
each iteration. The iterations converge to the solution for 
a given time step when the present and previous 
concentration values for each point in the entire domain 

k.i _ c J 



differ by less than an epsilon value (i.e., 
< £ for all points in the domain) . 



'b D 



'i/ D 



Figure (5) is a graphical representation of the program 
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Figure (5) - Formaldehyde Concentrations Within the 
Numerical Photo-reactor Model 
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output for selected values of time. The solution at t = 10 
sec is considered to be a steady state solution. This was 
determined by comparing the solutions at t = 10 sec and t = 
100 sec which did not differ to two decimal places for the 
vast majority of points in the domain. In addition, to 
ensure that a time-accurate, steady state solution was 
obtained, the computer program was run at At=0.005 s (half 
the original time step) out to a 10 sec solution. The 
results of this run did not differ from the original 10 sec 
solution to over 4 decimal places for the entire domain. 

The numerical concentration predictions for the "start- 
up" of the reactor make sense from a practical standpoint. 
The smallest concentrations are seen at the bulb surface and 
wake regions due to the low fluid velocities there. In 
addition, these regions of lesser concentration migrate 
further downstream with increased time as they are swept 
into the bulb wake region. 

The steady state isopleths also appear reasonable in 
form. As expected, these concentration contours generally 
form the same shape as the streamlines depicted in the 
FLUENT output due to the dominance of the convection terms. 
However, the isopleth shapes differ from the streamlines due 
to diffusive transport and photolysis. In the next chapter 
it will be shown by comparison with the simplified flowfield 
model that the magnitude of formaldehyde reductions shown in 
Figure (5) are reasonable predictions. 



Chapter 8 

ANALYSIS OF RESULTS 
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In order to obtain a general performance check on the 
numerical solution, the simplified flow model developed in 
Chapter 5 can be employed by inputting the numerical model 
parameters discussed in Chapter 7 into Equation [30] as 
follows : 

Co = 50 (PPM) 

K = 0.35 (1/s) 
v = 2 . 0 (m/ s) 

x = 0.2121 (m) 

Substituting values into equation [30] yields: 

C = Co exp{-Kx/v} 

= 50.0PPM exp{-0.35s _1 (0.2121m) /2.0m/s} 

= 48.17 PPM 



The above result compares favorably with the 49.37 PPM 
average steady state concentration leaving the numerical 
reactor model which was found through a mass average taken 
across the exit as follows: 



'exit, avg. 



ZviCi/Svi 



[45] 



In the next paragraph the analytical and numerical results 
are compared in terms of destruction efficiencies. 
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An additional comparison between the numerical and 
simplified flow field models is obtained by introducing 
parameters which define the destructive efficiency of the 
photolytic reactor. Since it is the objective of the 
reactor model to eliminate formaldehyde, the efficiency of 
the entire reactor section can be defined as follows: 



In addition to the above equation it is also useful to 
define an efficiency which is based on the bulb height 
alone. This term is defined as follows: 



Equation [47] calculates an efficiency based on a collimated 
region with a height equal to the radius of the bulb. In 
Equation [47] the term C e xit, avg. / c inlet is taken from the 
centerline of the reactor up to the bulb radius. Table-1 
compares the efficiency of the steady state numerical 
simulation of the photolytic reactor section with that of 
the simplified flow field model. 



E 



^exit, avg. ^inlet^ 



[46] 




[47] 



40 





Numerical Model 
(PPM) 


Simplified flow 
Field Model 
(PPM) 


C (average across 
entire reactor 
section) 


49.37 


48.17 


C (average up to 
bulb radius) 


47.37 


48.17 


E (based on entire 
reactor cross 
section) 


1.26% 


3.64% 


E B (based on bulb 
width) 


5.26% 


3.64% 



Table-1 Comparison of Numerical and Simplified Flow Models 



For this set of conditions it appears that the 
simplified flow field model is more efficient than the 
numerical reactor model when the efficiency is based on the 
entire reactor cross section. However, in examining the 
efficiencies based on the bulb cross section, the numerical 
reactor model which incorporates the actual flow conditions 
is superior to the simplified flow field model in its 
ability to destroy the formaldehyde. This result leads one 
to believe that a reactor containing a combination of bulbs 
which allows the incoming flow to "see" bulbs across the 
entire reactor area will out-perform the idealized columnar 
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reactor . 

The advantage of using multiple bulbs is best 
illustrated by considering a single bulb reactor of infinite 
height. In such a reactor the UV bulb would have little 
effect on the VOC concentrations in the far stream. As 
shown in Figure (5) , the bulb has the most influence 
(concentration reductions) in a columnar region which is one 
to two bulb diameters in height. In other words, the 
incoming flow which exists outside this columnar region 
tends to escape untreated. Therefore, an arrangement of 
bulbs which guarantees that all the incoming flow will 
encounter a bulb cross section (ie., a staggered stack of 
bulbs/columnar regions in parallel with no "escape" regions 
in between) would have a destructive efficiency of no less 
than that given by Equation [47] . 

In addition to the general advantages of a parallel 
bulb arrangement as described above, large-scale destructive 
efficiency improvements for both series and parallel 
arrangements are predicted by drawing an analogy between 
particle collection and VOC destruction in a photolytic 
reactor. For such an analogy the multi-bulb efficiency is 
given by the following (Heinsohn, 1991) : 

Overall = 1 " d ~ E B ) N [48] 



where, N is the number of bulbs. 
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A twenty-five (N = 25) bulb arrangement with single bulb 
efficiencies as given by Equation [47] (E B = 5.26% in this 
case) would have an overall destructive efficiency as 
follows : 

E overall = 1 “ U - E B ) N = 1 - (1 - 0.0526) 25 
= 74.1% 

Although the particle collector/photolytic reactor 
analogy drawn here predicts very optimistic results, it will 
be necessary to confirm this relationship with a multi-bulb 
numerical model and/or experimental testing to achieve a 
high degree of confidence. 



43 



Chapter 9 

CONCLUSIONS 



This analysis clearly demonstrates the influence of 
actual flow conditions on the effectiveness of a photolytic 
reactor. Despite the simplifying (but reasonable) 
assumptions made in development of the numerical model, it 
provides valuable insight into the inner workings of the 
actual photolytic reactor. The graphical results of Figure 
(5) show that the velocity gradients, recirculation and 
mixing which occur within the reactor contribute greatly to 
the overall reduction of reactants. Since recirculation and 
mixing greatly increase the reduction of formaldehyde in the 
reactor, one would expect that unsteadiness in the flow 
field would tend to improve the results (noting that 
unsteadiness exists in the actual flow field under 
consideration due to vortex shedding) . In addition, one 
would expect that an optimum Reynolds number, R e could be 
found (input velocity) which would maximize the benefits of 
mixing while still providing an acceptable level of air 
treatment. Therefore, future numerical optimization studies 
should include both an unsteady flow field analysis and a 
comparison of results for different Reynolds numbers. 

Although the simplified flow field model developed in 
this report can be used to obtain a rough estimate of 
reactor performance, it is unable to predict the effects of 
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the actual flow field which were found to be of great 
importance through an analysis of the numerical model 
results. This fact reinforces the need to use a numerical 
model in the optimization of a photolytic reactor design. 

Comparison of efficiencies between the simplified flow 
field and numerical models provides some highly instructive 
results. Surprisingly, though the light intensity in the 
numerical model is more than five times less than that of 
the simplified flow field model for a large portion of the 
domain (due to the radial light attenuation from the lamp 
surface) , the numerical model provides a comparable 
formaldehyde reduction performance when the efficiency is 
based on the entire reactor cross section. In addition, it 
appears that increasing the number of bulbs will greatly 
improve overall reactor effectiveness as indicated by a 
comparison of efficiencies based on the bulb width. Due to 
the optimistic results obtained in comparing the numerical 
and simplified flow field models, the next logical step for 
future research in this area should include the modeling of 
larger flow domains which contain a larger quantity of 
bulbs . 

This undertaking also showed that additional factors 
can readily be taken into consideration in future reactor 
modeling. The computer program of Appendix A and techniques 
described in this analysis can be extended to include 
additional chemical reactions and other factors such as 
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Beer-Lambert light attenuation by reactants which will 
result in a numerical model of much greater accuracy. 
Although this work is of limited scope for practical 
reasons, expansion of this analysis including the use of a 
high-speed supercomputer could achieve an accurate numerical 
simulation of the entire photolytic reactor of Figure (2) . 

Despite the good correlation between the numerical and 
simplified flow field models developed in this report, a 
high degree of confidence can only be achieved through 
laboratory experiments whose results can be compared to the 
numerical ones. Therefore, future work should include 
experimental verification. This will be especially 
important when additional reactants and other factors are 
introduced. In the future, the benefits of a closed form 
solution as experienced in this development (simplified flow 
field model) will not be possible when modeling a highly 
complex reactor. 
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COMPUTER PROGRAM 
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C SET UP VARIABLES 

IMPLICIT DOUBLE PRECISION (A-H,0-Z) 

DIMENSION VX (100, 33) ,VY(100,33) , FMC (10 0, 33) , FMCOLD (100 , 33) 

X0=4 .286 
Y0=0 . 0 
R0=1 . 50 
M=100 
N=31 

H=0. 002143 
EPSI=0. 0000001 
DT=0 . 01 
ITER-1 
ITERT=1 
ITEST=1 
I TERMAX=2 0 0 
I TERTMAX= 1 0 

C ENTER FORMALDEHYDE CONSTANTS 
FMRC0=0 . 35 
FMD=0 .00002 
FMC0=50 . 0 

C READ IN FLUENT GENERATED VELOCITIES 

OPEN (UNIT=5 , FILE=' SINGLEBULBLIS.DAT' , STATUS=' OLD' ) 

K=1 

KK=1 

KKK=10 

2 READ (5,5) ( (VX (1 , 33- J) , I=KK, KKK) , J=1 , N) 

K=K+1 

KKK=K*10 

KK=(K*10) -9 

IF (K.LE.10) GO TO 2 

K=1 

KK=1 

KKK=10 

3 READ (5,5) ( (VY (I , 33- J) , I=KK, KKK) , J=1 , N) 

K=K+1 

KKK=K*10 
KK= (K*10 ) -9 
IF (K.LE.10) GO TO 3 
CLOSE (UNIT-5) 

5 FORMAT (5X, 10E12 .4) 

C SET THE INITIAL FORMALDEHYDE CONCENTRATIONS 
DO 10 1=1, M 
DO 10 J=1 , N+2 
FMC (I, J) =FMC0 
FMCOLD ( I, J) =FMC0 
10 CONTINUE 

C CALCULATE THE CONCENTRATIONS AT EACH POINT USING GAUSS- SEIDEL METHOD 
15 DO 50 1=2, M-l 
DO 50 J=2 , N+l 
TEMP=FMC ( I , J) 

C CALCULATE THE PHOTOLYTIC RATE CONSTANT 

FMRC=FMRC0 *R0 / ( (100. *H* (1-1) -X0) **2 . + (100 . *H* (J-2) ) **2 .+ 

$0 . 00000000001) **0.5 
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C APPLY BOUNDARY CONDITIONS (BC'S) AT THE BULB (J=2, 1=14; J=4, 1=15; 

C J=6 , 1=16, J=7, 1=17, J=8, 1=19) 

IF (( (J.EQ.2) .AND. (I.EQ.14) ) .OR. ( (J.EQ.4) .AND. (I.EQ.15) ) .OR. 

$ ( (J.EQ.6) .AND. (I.EQ.16) ) .OR. ( (J.EQ.7) .AND. (I.EQ.17) ) .OR. 

$ ( (J.EQ.8) .AND. (I.EQ.19) ) ) 

$FMC (I, J) = (1/ (4 . *FMD/H**2 . +FMRC+1 . /DT) ) * 

$ ( (FMD/H**2 . ) * (2.*FMC(I-1, J)+2.*FMC(I, J+l) ) +FMCOLD ( I , J) /DT) 

C APPLY BOUNDARY CONDITIONS (BC'S) AT THE BULB ( J=2 , 1=28 ; J=4 , 1=27 ; 

C J=6, 1=26, J=7 , 1=25, J=8, 1=23) 

IF (( (J.EQ.2) .AND. (I .EQ.28) ) .OR. ((J.EQ.4) .AND. (I .EQ.27) ) .OR. 

$ ( (J.EQ.6) .AND. (I.EQ.26) ) .OR. ( (J.EQ.7) .AND. (I.EQ.25) ) .OR. 

$ ( (J.EQ.8) .AND. (I.EQ.23) ) ) 

$FMC (I, J) = (1/ (4 . *FMD/H**2 . +FMRC+1 . /DT) ) * 

$ ( (FMD/H**2 . ) * (2 . *FMC (1+1, J) +2 . *FMC(I, J+l) ) +FMCOLD (I, J) /DT) 

C SKIP AREAS INSIDE THE BULB 

IF ( ( ( (J.GE.2) .AND. (J.LE.3) ) .AND. ( (I.GE.15) .AND. (I.LE.27) ) ) .OR. 

$ ( ( (J.GE.4) .AND. (J.LE.5) ) .AND. ( (I.GE.16) .AND. (I.LE.2 6) ) ) .OR. 

$ ( (J.EQ. 6) .AND. ( (I.GE.17) .AND . (I.LE.25) ) ) .OR. ( (J.EQ.7) .AND. 

$ ( (I.GE.18) .AND. (I.LE.24) ) ) .OR. ( (J.EQ.8) .AND. ( (I.GE.20) .AND. 

$ (I.LE.22) ) ) ) GO TO 50 
C CALCULATION FOR GENERAL POINTS 
IF ( (J.GE. 9) .OR. 

$ ( ( (I.LE. 13) .OR. (I.GE.29) ) .AND. (J.EQ.2) ) .OR. 

$ ( ( (I.LE.14) .OR. (I.GE.28) ) .AND. ( (J.GE. 3) .AND. (J.LE.4) ) ) .OR. 

$ (( (I.LE. 15) .OR. (I.GE.27) ) .AND. ( (J.GE. 5) .AND. (J.LE.6) ) ) .OR. 

$ (( (I.LE. 16) .OR. (I.GE.26) ) .AND. (J.EQ.7) ) .OR. 

$ (( (I.LE. 18) .OR. (I.GE.24) ) .AND. (J.EQ.8) ) ) 

$FMC (I, J)= (1/ ( (ABS (Vx (I, J) ) +ABS (VY (I, J) ) +4 . *FMD/H) /H+FMRC+1 . /DT) ) * 

$ ( ( (ABS (VX(I, J) ) +VX(I, J) ) *FMC (1-1, J) + (ABS (VX (I, J) ) -VX (I, J) ) * 

$FMC (1 + 1, J) + (ABS (VY(I, J) )+VY(I, J) ) *FMC (I, J-l) + 

$ (ABS (VY (I, J) ) -VY (I, J) ) *FMC(I, J+l) ) / (2.*H) + (FMD/H**2. ) * 

$ (FMC(I-1, J)+FMC(I+1, J)+FMC(I, J-1)+FMC(I, J+l) ) +FMCOLD (I, J) /DT) 

IF (ABS (FMC (I, J) -TEMP) .GT.EPSI) ITEST=ITEST*0 . 

50 CONTINUE 

C ASSIGN VALUES REQUIRED DUE TO BOUNDARY CONDITIONS 
DO 52 1=2, M-l 
FMC (1,1) =FMC (1,3) 

FMC (1,33) -FMC (I, 31) 

52 CONTINUE 

FMC (20, 8 ) =FMC (20, 10) 

FMC (21, 8) -FMC (21, 10) 

FMC (22, 8) =FMC (22, 10) 

FMC (18, 7) =FMC (18, 9) 

FMC (24, 7 ) =FMC (24, 9) 

FMC (16, 5) -FMC (14, 5) 

FMC (26,5) =FMC (28,5) 

FMC (15,3) =FMC (13,3) 

FMC (27,3) =FMC (29,3) 

C APPLY BC'S AT END (DC/DX=0) 

DO 55 J=1 , N+2 
55 FMC (M, J) =FMC (M-l, J) 

IF (ITER.GE. ITERMAX) GO TO 60 
IF (ITEST.EQ. 1) GO TO 60 
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ITER=ITER+1 
ITEST=1 
GO TO 15 

60 IF (ITERT.EQ. ITERTMAX) GO TO 64 
DO 62 1=2, M 
DO 62 J=2 , N+l 
62 FMCOLD ( I , J) =FMC ( I , J) 

ITER=1 

ITERT=ITERT+1 
ITEST=1 
GO TO 15 

C WRITE OUT RESULTS 

64 OPEN (UNIT=6 , FILE=' NEWSINGLE . OUT ' , STATUS='NEW' ) 

OPEN (UNIT=8,FILE=' SINGLEPLOT. DAT' , STATUS=' NEW' ) 

DO 100 1=15,27 

DO 100 J=2, 8 

C ASSIGN VALUES TO ACHIEVE A DECENT CONTOUR PLOT USING GNUPLOT 

100 IF ( ( ( (J.GE.2) .AND. (J.LE.3) ) .AND. ( (I.GE.15) .AND. (I.LE.27) ) ) .OR. 
$ ( ( (J.GE.4) .AND. (J.LE.5) ) .AND. ( (I.GE.16) .AND. (I.LE.26) ) ) .OR. 

$ ( (J.EQ.6) .AND. ( (I.GE.17) .AND. (I.LE.25) ) ) .OR. ( (J.EQ.7) .AND. 

$ ( (I .GE. 18) .AND. (I.LE.24) ) ) .OR. ( (J.EQ. 8) .AND. ( (I.GE.20) .AND. 

$ (I.LE.22) ) ) ) FMC (I, J) =30 . 

SUM=0 . 

C CALCULATE EXIT AVERAGES 
SUMD=0 . 

SUM1=0 . 

SUM1D=0 . 

DO 80 J=2, N+l 
SUMD=SUMD+VX (100, J) 

80 SUM=SUM+FMC (100, J) *VX(100, J) 

SUM=SUM/SUMD 

DO 81 J=2 , 9 
SUM1D=SUM1D+VX ( 1 00 , J) 

81 SUM1=SUM1+FMC (100, J) *VX (100, J) 

SUM1=SUM1/SUM1D 

K=1 

KK=1 

KKK=10 

WRITE (6, 71) ITER, ITERT, SUM, SUMl,DT,EPSI 

65 WRITE (6, 70) ( (FMC (I, 33-J) , I=KK,KKK) , J=1,N) 

K=K+1 

KKK=K*10 
KK= (K*10) -9 
IF (K.LE. 10) GO TO 65 
DO 67, J=2,N+1 
DO 66, 1=1, M 

66 WRITE (8, 68) FMC(I,J) 

WRITE (8, 69) 

67 CONTINUE 

68 FORMAT (F12 . 8) 

69 FORMAT (IX) 

70 FORMAT (10F12 . 8) 

71 FORMAT (216, 4E12. 4) 

STOP 

END 



Appendix B 

ATMOSPHERIC CHEMISTRY 
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WELL-MIXED MODEL OF 

PHOTO-OXIDATION AND RADICAL OXIDATION PROCESSES 



Contents: 



Photolysis Rates for 337 nm (laser) 

Photolysis Rates for 253 & 185 nm. (UV bulb) 

Photolysis Rates for laser & UV bulb (combo) 



Radical Chemistry 

Reactions: 

Formaldehyde 

Methanol 
Carbon Dioxide 

Water Vapor 

Ozone 

Atomic Oxygen 

NO 

H2O2 

N02 

HONO 

Hydroxyl Radicals..(OH) 

Hydro peroxyl Radicals.. (HO2) 
Carbon Monoxide (CO) 

Atomic Hydrogen..(H) 

Formyl Radical (HCO) 

Atmospheric Gases 

References. 



2 

3 

4 

5 



6 



/ 



8 

9 

10 

11 

12 

13 

14 

15 



Note 

R represents a reaction; P represents a photolvtic reaction; KR 
represents a rate constant; KP represents a photolytic rate constant. 
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Photolysis: 

X for the laser experiment will be 337 nm 
input file "RATE" 

Ephoton= 5 90* 1 O' 2 1 J/Photon 

Actinic Flux=20*10 15 Photons/(cm 2 *s) 



HCHO+hv -* H*+HCO* 


X<300nm 


4.08 E-5 s* 1 


fPlAJ 


Atkinson 


HCHO+hv — COH 2 


Xa300nm 


2.97E-4 s' 1 


[PI B] 


Atkinson 


CH 3 OH+I 1 V — • negligible 




0 


[P2] 




Q3+hv-*0(*D)402 


X<;41 lnm 


0 


[P3A] 


Atkinson 


03 +hv — 0(3p)+0 2 


X<1180nm 


0, (<J> = 0 ) 


[P3B] 


Atkinson 


02+hv -*0( 1 D)+0(3p) 


X<175nm 


0 


[P4AJ 


Atkinson 


02+hv -* 0(3p)+0(3p) 


X<242nm 


0 


[P4B] 


Johnson/ 

Atkinson 


02+hv-0(lD)+0(lD) 


X<137nm 


0 


[P4C] 


Johnson/ 

Atkinson 


H 2 02+hv-*0H*+0H* 


X<557nm 


1.40E-5 s-1 


[PSA] 


Atkinson 


H 2 02+hv — H 2 0+0( l D) 


X<359nm 


0, (<t> = 0 ) 


[P5B] 


Atkinson 


H 2 0 2 +hv -H*+H02* 


X<324nm 


0, (<J> = 0) 


[P5C] 


Atkinson 


H202+hv -*2H*+02 


X<213nm 


0 


[PSD] 


Atkinson 


N02+hv-N0+O(lD) 


Xs244nm 


0 


[P 6 A] 


Atkinson 


NCb+hv — N0+0( 3 P) 


X<398nm 


7.08E-3 s' 1 


[P 6 B] 


Atkinson ' 


H 2 0+hv -H 2 +0( 3 P) 


X<243nm 


0 


f p 7A] 


Atkinson 


H 2 0+hv —H*+OH* 


X<239nm 


0 


[P/B] 


Atkinson 


H 2 0+hv — H 2 +0( 1 D) 


X<176nm 


0 


[P7C] 


Atkinson 


C0 2 +hv — CO+0( *D) 




negligible 






COz+hv — C0+0( 3 P) 




negligible 






HONO+hv — NO+OH* 


X<591nm 


3.23E-3S-1 


[P50A] 


Atkinson 


HONO+hv — H*+NOz 


X<367nm 


0 , (4> = 0 ) 


[P50B] 


Atkinson 


HONO+hv -*HN0*+0( 3 P) 


X<283nm 


0 


[P50C] 


Atkinson 



Photolysis: 

X for the UV Bulb experiment will be 253 nm and 185 nm 
input file ’•BULB” 

Hphoton=5 90* 10-21 J/Photon 

Actinic Flux: J(253)=1.14*10 16 Photons/(cm 2 *s) 
J(185)=3.52*lQl 4 Photons/(cm 2 *s) 



HCHCM-hv-H*+HCO* 


X<300nm 


1.30E-5s - 1 


[P1A] 


Atkinson 


HCHO+hv — CCh-H 2 


X^300nm 


2.16E-5 s ' 1 


[PI BJ 


Atkinson 


CH 3 OH+I 1 V — negligible 






[P2] 




03+hv-0(iD)+0 2 


Xs411nm 


1.12E-1 s ' 1 


[P3A] 


Atkinson 


03 +hv -*• 0( 3 P)+0 2 


X<1180nm 


1.32E-2 s’l - 


- [P3B] 


Atkinson 


02+hv^0(!D)+0(3p) 


X<175nm 


0 


[P4A] 


Atkinson 


0 2 +hv — 0( 3 P)+0( 3 P) 


X<242nm 


1.23E-9 s*l 


[P4B] 


Johnson/ 

Atkinson 


02 +hv- 0 (lD)+ 0 (lD) 


X<137nm 


0 


[P4C] 


Johnson/ 

Atkinson 


H 2 0 2 +hv -»> OH*+OH* 


X<557nm 


8.38E-4 s ' 1 


[PSA] 


Atkinson 


H 2 C> 2 +hv — H 2 0+0(!D) 


X<359nm 


0 ,(<J> = 0 ) 


[P5B] 


Atkinson 


H 2 02 +hv — *-H*+H0 2 * 


X<324nm 


0 , (<t> = 0 ) 


[P5C] 


Atkinson 


H 2 02 +hv -»• 2H*+0 2 


X<213nm 


0 , (<P = 0 ) 


[PSD] 


Atkinson 


NCb+hv -»• N0+0(lD) 


Xs244nm 


0 , (O = 0 ) 


[PGA] 


Atkinson 


N 02 +hv^N 0 + 0 ( 3 p) 


X<398nm 


2.47E-4 s * 1 


[P 6 B] 


Atkinson 


H 2 0 +hv — ►H 2 +0( 3 P) 


X<243nm 


0 ,(O = 0 ) 


[P7A] 


Atkinson 


H 2 0+hv — H*+OH* 


X<239nm 


1.09E-5 s"l 


[P7B] 


Atkinson 


H 2 0+hv — H 2 +0(iD) 


X<176nm 


0 


[P7q 


Atkinson 


C 02 +hv^COO(iD) 




negligible 






C02+hv^C0+0( 3 P) 




negligible 






HONO+hv — NO+OH* 


X<591nm 


1.69E-3 s * 1 


[P50A] 


Atkinson 


HONO+hv — H*+N0 2 


X<367nm 


0 , (O = 0 ) 


[PS OB] 


Atkinson 


HONO+hv -*HN0*+0( 3 P) 


X<283nm 


0 ,(O = 0 ) 


[P50C] 


Atkinson 



Photolysis: 
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X for the UV Bulb and the laser experiment will be 337, 253, and 185 nm 
input file "COMBO" 

Ephoton=590*10-2l J/Photon 

Actinic Flux: J(337)=20.0*10 15 Photons/(cm 2 *s) 

J(253)=1.14*10 16 Photons/(cm 2 *s) 

J(185)=3.52*10 14 Photons/(cm 2 *s) 



HCHO+hv -* H*+HCO* 


X<300nm 


1.30E-5S - 1 


[P1A] 


Atkinson 


HCHO+hv -* G>H 2 


>«i300nm 


3.86E-5 s ' 1 


(P1B] 


Atkinson 


CH 30 H+hv — negligible 






[P2] 




Q 3 +hv- 0 (iD)+ 0 2 


X^41 lnm 


1.12E-1 s * 1 


(P3A] 


Atkinson 


O 3 +hv — ► 0(3 P) +0 2 


X<1180nm 


1.32E-2 s * 1 


[P3B] 


Atkinson 


02 +hv -* 0( 1 D)+0(3p) 


X<175nm 


0 


[P4A] 


Atkinson 


02 +hv -* 0(3p)+0(3p) 


X<242nm 


1.23E-9 s * 1 


[P4B] 


Johnson/ 

Atkinson 


02+hv -0(lD)+0( ID) 


X<137nm 


0 


[P4C] 


Johnson/ 

Atkinson 


H 2 0 2 +hv — OH*+OH* 


X<557nm 


8.52E-4 s ' 1 


[P5AJ 


Atkinson 


H 2 02+hv — H 2 G+0( 1 D) 


X<359nm 


0 , (0 = 0 ) 


[P5B] 


Atkinson 


H 2 02 +hv -*H*+H02* 


X<324nm 


0 , (O = 0) 


[P5C] 


Atkinson 


H 2 02+hv -*• 2H*+0 2 ' 


X<213nm 


0 , (0 = 0) 


[PSD] 


Atkinson 


N02+hv -* N0+0( 1 D) 


Xs244nm 


0 , (0 = 0) 


[P6A] 


Atkinson 


N 02 +hv — N0+0(3p) 


X<398nm 


7.42E-4 s * 1 


[P6B] 


Atkinson 


H 2 0+hv -H 2 +0(3p) 


X<243nm 


O ,(0 = O) 


[P7A] 


Atkinson 


H 2 0+hv — H*+OH* 


X<239nm 


1.09E-5 s ' 1 


[P7BJ 


Atkinson 


H20+hv— H 2 +0(1D) 


X<176nm 


0 


[P7q 


Atkinson 


C02+hv— C0+0(1D) 




negligible 






C02+hv -* C0+0(3p) 




negligible 






HONO+hv — NO+OH* 


X<591nm 


2.6 IE- 3 s ’ 1 


[P50A] 


Atkinson 


HONO+hv — H*+N0 2 


>.<367nm 


0, (0 = 0) 


[P50B] 


Atkinson 


HONO+hv -*HN0*+0(3p) 


X<283nm 


0, (0 = 0) 


[P30C] 


Atkinson 
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Radical Chemistry: 



0H*+Q3=H02*+02 


6.7*10" 14 


cm 3 molecule"! s'! 


[Rl] 


DeMore 


H02*+O3=0H*+202 


2.0*10" J5 


cm 3 molecule"! s'! 


(R2] 


DeMore 


NO+O3=N02+O2 


1.8*10" 14 


cm 3 molecule* 1 s"! 


[R3] 


DeMore 


N02+03=NQ34O2 


1.2*10" 13 


cm 3 molecule"! s'! 


[R4] 


DeMore 


HC0*4O2=H02*4<X) 


5.6*10* 12 


cm 3 molecule"! s"! 


[R9] 


Anatasi 


HCHO+OH*=H*4CO+H20 


l.l*10" n 


cm 3 molecule"! s"! 


[R10] 


Seinfeld 


0H*+OH*=H 2 O+O(iD) 


1.9*10" 12 


cm 3 molecule"! s"! 


[Rill 


DeMore 


0H*+H02*=H 2 O+O2 


1.1*10" 10 


cm 3 molecule"! s"! 


[R12] 


Atkinson 


0H*+H 2 02=H2O+H02* 


1.7*10*12 


cm 3 molecule"! s"! 


[R13] 


DeMore 


N02+0H*=IiN03 


l.ino- 1 1 


cm 3 molecule"! s"! 


[R14] 


Seinfeld 


C0+0H*=C02+H* 


2.9*10*! 3 


cm 3 molecule"! s"! 


[R1S] 


DeMore 


N+OH*=NO+H* 


4.9*10*! 1 


cm 3 molecule"! s'! 


[R16] 


Atkinson 


hcho+oh*=h 2 ohgo 


1.0*10*!! 


cm 3 molecule"! s"! 


[R17] 


Atkinson 


H*+0? = HOz* 


1.2*10*1 2 


cm 3 molecule" 1 s'! 


[Rl 83 


Anatasi 


H*+03 = OH* + 02 


2.9*10*11 


cm 3 molecule"! s"! 


[R19] 


Anatasi 


H02*+H02*=H 2 0 2 +0 2 


1.6*10-12 


cm 3 molecule"! s" 1 


[R20] 


DeMore 


NO+H02*=N02+OH* 


8.3*10*12 


cm 3 molecule"! s"! 


[R22] 


DeMore 


N02+H0 2 *=H0N0+O2 


1.4*10*1 4 


cm 3 molecule"! s"! 


[R23] 


Stockwell 


CC^H02*=C0 2 4OH* 


1.9*10* 32 


cm 3 molecule"! s"! 


[R24] 


Atkinson 


0(1D)+N 2 =0( 3 P)+N 2 


1.8*10-! 1 


cm 3 molecule"! s*! 


[R30] 


DeMore 


0( 3 P)+O2+M=O 3 4M 


2.9*10*11 


cm 3 molecule"! s"! 


[R31] 


Seinfeld 


NO+O( 3 P)+M=N02-fM 


3.0*10-11 


cm 3 molecule"! s"! 


[R32] 


Atkinson 


NO2+O( 3 P)=O2+N0 


9.7*10*12 


cm 3 molecule"! s"! 


[R33] 


DeMore 


0( 3 P)+0 3 =02407 


9.5*10*1 5 


cm 3 molecule"! s"! 


[R35] 


DeMore 


0( 1 D)+H2O=0H*+0H* 


2.2*10*1° 


cm 3 molecule" 1 s'! 


[R36] 


DeMore 


0( lD)+O 3 =02+0( 3 P)+0( 3 P) 


1.2*10-H 


cm 3 molecule"! s"! 


[R37] 


DeMore 


0(lD)+02=0(" 3 P)+02 


4.0*10*1 1 


cm 3 molecule" 1 s"! 


[R38] 


DeMore 


N+0 2 =NO+O(3p) 


8. 9*10-1 7 


cm 3 molecule" 1 s" 1 


[R39] 


Atkinson 


N+0^=N0+O2 


1 . 0 * 10 - 1 ° 


cm 3 molecule"! s" 1 


[R40] 


Atkinson 


NO+OH*=HONO 


1 . 1 * 10-11 


cm 3 molecule"! s"! 


[R51] 


Atkinson 


H0N0+0H*=H 2 O+-NCb 


4.9*10*12 


cm 3 molecule"! s"! 


[R52] 


Atkinson 


CH 3 OH+OH*=0.85 |CH 2 OH*]+ 


7.8*10*1 3 


cm 3 molecule"! s" 1 


[R60] 


Atkinson 


0.15|CH 3 O*-rH2OJ 










CH 2 0H*+0?=HCHOfH0 2 * 


9.8*10*12 


cm 3 molecule"! s" 1 


[R61] 


Atkinson 


CH 3 0*+02=HCH04-H02* 


1.9*10*1 5 


cm 3 molecule" 1 s* 1 


[R62] 


Atkinson 



Formaldehyde Reactions (HCHO): 
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Photolysis: 



HCHO+hv —H*+HCO* 


X<300nm 


1.30E*5s *1 


[PI A] 


HCHO+hv — CO+H 2 


Xi300nm 


3.86E-5 s * 1 


[P1B] 


Radical Chemistry: 


HCHO+OH*=H*+O>H20 


l.l* 10 * n cm 3 


molecule*! s*! 


[RIO] 


HCHC4OH*=H 2 OfHC0* 


1 . 0 * 10 * n cm 3 


molecule*! s*! 


[R17] 


CH 2 0H*+O2=HCHO+H02* 


9.8*10 * 12 cm 3 


molecule ' 1 s * 1 


[R61] 


CH 3 0*+02=HCHO+H02* 


1.9*10 * 15 cm 3 


molecule*! s*! 


[R62] 



d[HCHOJ/dt=-PlA-PlB-R10-R17+R61+R62 



Methanol Reactions (CH 3 OH): 

Photolysis: 

CH 30 H-rhv — negligible [P2] 

Radical Chemistry: 

CH 3 OH+OH*=0.85 {CH 2 OH*l+ 7.8*10*13 cm 3 molecule*! s * 1 [R60] 

0.15{CH 3 O*+H2O| 

d[CH 3 OH]/dt=-R60 



Carbon Dioxide Reactions (CO 2 ): 

Radical Chemistrv: 

CO + OH*-** CO 2 + H* 2.2*10" 13 cm3 molecule*! s*! [R15] 

CO + HO 2 *-* CO 2 + OH* 1.9*10*32 cm3 molecule*! s*! [R24] 



Seinfeld 

Seinfeld 



Seinfeld 

Atkinson 

Atkinson 

Atkinson 



Atkinson 



DeMore 

Atkinson 



d[C02]/dt=P15+P24 



Atomic Oxygen Reactions 


(0): 








61 


Photolvsis: 












03+hv-O(iD)4O2 


X*41 lnm 


1.12E-1 s ' 1 


[P3A] 


Atkinson 


03 +hv -*0(3p)+02 


X<1180nm 1.32E-2S- 1 


[P3B] 


Atkinson 


02+hv -*0( 1 D)+0(3p) 


X<175nm 


0 


[P4A] 


Atkinson 


02 +hv- 0 ( 3 p)+ 0 ( 3 p) 


X<242nm 


1.23E-9 s ' 1 


[P4BJ 


Atkinson 


02 +hv -O(lD)+0(lD) 


X<137nm 


0 


[P4C] 


Atkinson 


N02+hv— NG+0(1D) 


Xs244nm 


0, (0 = 0) 


[PGA] 


Atkinson 


N 02 +hv-N 0 + 0 ( 3 p) 


X<398nm 


7.42E-4 s ' 1 


[PGB] 


Atkinson 


H 2 0+hv -H 2 +0(3p) 


X<243nm 


0, (0 = 0) 


[P7A] 


Atkinson 


H 2 0+hv-H2+0( 1 D) 


X<176nm 


0 


[P7C] 


Atkinson 


H 2 02+hv -* H 2 0+0( 1 D) 


X<359nm 


0, (0 = 0) 


[P5B] 


Atkinson 


HONO+hv -HN0*+0(3P) 


X<283nm 


0, (0 = 0) 


[P50C] 


Atkinson 


Radical Chemistrv: 
0H*+0H*=H2O+O(iD) 


1.9*10 ' 12 


cm3 


molecule ' 1 s ' 1 a 


[Rll] 


DeMore 


0(1D)+N 2 =0(3P)+N2 


l. 8 * 10- ]1 


cm3 


molecule ' 1 s ' 1 


[R30] 


DeMore 


0(3P)+02+M=0 3 +M 


2.9*10 ' 1 1 


cm3 


molecule ' 1 s ' 1 


[R31] 


Seinfeld 


NO+0(3P)+M=N02+M 


3.CT10 ' 11 


cm3 


molecule ' 1 s ' 1 


[R32] 


Atkinson 


N02+0( 3 P)=02+n6 


9.7*10 ' 12 


cm3 


molecule ' 1 s ' 1 


[R33] 


DeMore 


0(3p)+0 3 =20 2 


9.5*10 ' 15 


cm3 


molecule ' 1 s ' 1 


[R35] 


DeMore 


0 ( 1 D)+H 2 O=20H* 


2 . 2 * 10' 10 


cm3 


molecule ' 1 s ' 1 


[R36] 


DeMore 


0(1D)+O 3 =02+20(3P) 


1 . 2 *lO' n 


cm3 


molecule ' 1 s * 1 


[R37] 


DeMore 


0( 1 D)+02=0(3p)+02 


4.0M0 ' 11 


cm3 


molecule ' 1 s ' 1 


[R38] 


DeMore 


N+02=N0+0(3p) 


8.9*10 - 17 


cm3 


molecule ' 1 s ' 1 


[R39] 


Atkinson 



a Assumed the products of Rll were negligible. 
Use Steady State Assumption: . 



[O(3p)] ss =(P3B+P4A+P4B+P6B+P7A+P50C+R30+R37+R38+R39)/(KR33* 
N02*KR35*03+KR32*N0+KR3 1*02) 

[0( 1 D)] SS =(P3A+P4A+P4C+P6A+P7C+P5B+R1 1 )/(KR30*N2+KR36*H2O+ 
KR37*03+KR38*02) 



NO Reactions: 
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Photolysis: 



N02+hv-NO+O(lD) 


Xs244nm 


0, (<i> = 0) 


[P6A] 


Atkinson 


N02+hv-N0+0(3p) 


X<398nm 


7.42E-4 s' 1 


[P6B] 


Atkinson 


HONO+hv — NO+OH* 


X<591nm 


2.61E-3 s' 1 


[P50A] 


Atkinson 


Radical Chemistrv: 


N0 + 03 = N02 + 0 2 


1.8*10' 14 cm 3 


molecule' 1 s' 1 


[R3] 


DeMore 


N + OH* = NO + H* 


4.9*10' n cm 3 


molecule' 1 s' 1 


[R16] 


Atkinson 


NO + H0 2 * = NO 2 + OH* 


8.3*10' 12 cm 3 


molecule' 1 s' 1 


[R22] 


DeMore 


•NO + 0(3P) + M=N02 + M 


3.0*10 _1 1 cm 3 


molecule' 1 s' 1 


[R32] 


Atkinson 


NO 2 + 0( 3 P) =02 + NO 


9.7*10' 12 cm 3 


molecule' 1 s' 1 


[R33] 


DeMore 


N+02=N0+O(3 P) 


8.9*10' 17 cm 3 


molecule' 1 s' 1 


[R39] 


Atkinson 


N+Q3=N0+02 


1.0*10' 16 cm 3 


molecule' 1 s' 1 


[R40] 


Atkinson 


NO+OH*=HONO 


l.mo- 11 cm 3 


molecule' 1 s' 1 


[R51] 


Atkinson 



d [NO]/ dt=P6A+P6B+P5 0A-R3 +R1 6-R22-R3 2+R3 3 +R3 9+R40-R5 1 



Hydrogen Peroxide Reactions (H2O2): 



Photolysis: 



H 2 O 2 -r hv— OH* + OH* 


X<577 nm 


8-52H-4 s’ 1 


[PSA] 


Atkinson 


H 2 O 2 + hv— H 2 0 + 0(!D) 


>.<359 nm 


0, (0 = 0) 


[P5B] 


Atkinson 


H 2 O 2 + hv— H*+H0 2 * 


X<324 nm 


0, (0 = 0) 


[PSC] 


Atkinson 


H 2 O 2 + hv— 2H*+02 


X<213 nm 


0, (0-0) 


[PSD] 


Atkinson 


Radical Chemistry: 


0H* + H 2 02 = H?0 + H0?* 


1.7*10' 12 cm 3 


molecule' 1 s' 1 


[R13) 


DeMore 


H02* + H02*=H 2 0 2 +02 


1.6*10' 12 cm 3 


molecule' 1 s' 1 


[R20] 


DeMore 



d [H 2 O 2 ] / d t=-P5 A-P5 B-P5 C-P5 D-R 1 3 +R20 



NO 2 Reactions: 
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Photolysis: 



N02+hv-*N0+0(lD) 


Xs244nm 


0, (0 = 0) 


[PGA] 


Atkinson 


NC>2+hv^NO+0(3p) 


X<398nm 


7.42E-4 s* 1 


[P6B] 


Atkinson 


HONO+hv— H*+N02 


k<367nm 


0, (O = 0) 


[P50B] 


Atkinson 


Radical Chemistry: 
NO+O3-N07+02 


1.8*10‘ 14 


cm 3 


molecule*! s*! 


[R3] 


DeMore 


N02403=N03+02 


1.2*10* 13 


cm 3 


molecule*! s'! 


m 


DeMore 


NO2+OH*=HN0 3 


l.l*10* n 


cm 3 


molecule* 1 s*! 


[R14] 


Seinfeld 


NO+-H0 2 *=N02+0H* 


8.3*10*! 2 


cm 3 


molecule*! s*! 


[R22] 


DeMore 


N02+H02*=H0NCMO2 


1.4*10* 12 


cm 3 


molecule*! s*! 


[R23] 


Stockwell 


NO+O( 5 P)+M=N02+M 


3-ono* 1 1 


cm 3 


molecule* 1 s* 1 


[R32] 


Atkinson 


N02+O(3p)=02+N0 


9.7*10* 12 


cm 3 


molecule*! s'! 


[R33] 


DeMore 


H0N0+0H*=H 2 0+NC)2 


4.9*10-12 


cm 3 


molecule*! s"! 


[R52] 


Atkinson 



d[N02]/dt=-P6A-P6B+P5 0B+R3-R4-R1 4+R2 2-R23 +R3 2-R3 3 +R5 2 



Nitrous Acid Reactions (HONO): 



Photolysis: 



HONO+hv —NO+OH* ' 


X<591nm 


2.6 IE- 3 s*l 


[P50A] 


Atkinson 


HONO+hv— H*+N02 


>.<3G7nm 


0, (4> = 0) 


[P50B] 


Atkinson 


HONO+hv-HNO*+C( 3 P) 


>.<283nm 


0 , (O = 0) 


[P50C] 


Atkinson 


Radical Chemistry: 


N0 2 +H02*=H0NO02 


1.4*10* 14 cm 3 


molecule*! s*! 


[R23] 


Stockwell 


NO+OH*=HONO 


l.l*10* n cm 3 


molecule*! s’! 


[R51] 


Atkinson 


HQN040H*=H 2 0+N02 


4.9*10*1 2 cm 3 


molecule*! s*! 


[R52] 


Atkinson 



d[HONO]/dt=-P50A-P50B-P50C+R23+R51-R52 



Hydroxyl Radical Reactions (OH*): 



Photolysis: 



H202+hv-0H*+0H* 


X<557nm 8.52E-4 s* 1 


[P5A] 


Atkinson 


H 2 0+hv— H*+OH* 


X<239nm 1.09E-5 s’ 1 


[P7B] 


Atkinson 


HONO+hv — NO+OH* 


X<591nm 2.61E-3 s* 1 


[P50A] 


Atkinson 


Radical Chemistry: 










0H*4O3=H0 2 *+O 2 


6.7*1044 


cm 3 molecule*! s' 1 


[Rl] 


DeMore 


H02*+03=0H*+202 


2.0*1045 


cm 3 molecule'! s"! 


[R2] 


DeMore 


HCH0+0H*=H*+C0+H 2 0 


1.1*1041 


cm 3 molecule*! s*! 


[R10] 


Seinfeld 


0H*+OH*=H 2 O+O( 1 D) 


1.9*1042 


cm 3 molecule*! s*! b 


[Rll] 


DeMore 


0H*+H0 2 *=H20+O2 


1.1*1040 


cm 3 molecule*! s'! 


[R12] 


Atkinson 


0H*+H 2 02=H20+H02* 


1.7*104 2 


cm 3 molecule*! s' 1 


[R13] 


DeMore 


N02+OH*=HN03 


1.1*1041 


cm 3 molecule*! s*! 


[R14] 


Seinfeld 


CO+0H*=CO2+H* 


2.2*104 3 


cm 3 molecule'! s' 1 


[R15] 


DeMore 


N+OH*=NO+H* 


4.9*104 1 


cm 3 molecule'! s'! 


[Rl<3] 


Atkinson 


H*+03 = 0H*+02 


2.9*1041 


cm 3 molecule*! s*! 


[R19] 


Anatasi 


N0+H02*=N02+0H* 


8.3*104 2 


cm 3 molecule*! s*! 


[R22] 


DeMore 


COH02~*=C02+0H* 


1.9*10* 32 


cm 3 molecule*! s*! 


[R24] 


Atkinson 


0(lD)+H 2 O=0H*+OH* 


2.2*1040 


cm 3 molecule*! s'! 


[R36] 


DeMore 


NO+OH*=HONO 


1.1*1042 


cm 3 molecule* 1 s* 1 


[R51] 


Atkinson 


H0N04OH*=H 2 O4-N02 


6.6*104 2 


cm 3 molecule'! s'! 


[R52] 


Atkinson 


CH 3 OH+OH*=0.85 fCH 2 OH*}+ 


7.8*1043 


cm 3 molecule'! s*! 


[R60] 


Atkinson 



0.15{CH 3 O*+H2O{ 



b Assumed [OH*] to be low, considered Rll to be negligible. 

Use Steady State Assumption: 

[OH*]ss=(P5A+P7B+P50A+R2+R19+R22+R24+2’ r R36)/(KRl*O3-i-KR10* 
HCH0+KR12*H02+KR13*H202+KR14*N02+KR15*C0+KR16*N+ 
KJR5 l*NO+KR52*HONO+KR60*CH3OH) 



Hydroperoxyl Radical Reactions (H0 2 *): 
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Photolysis: 



H 2 0 2 +hv-*H*+H0 2 * 


X<324nm 0, (4>=0) 


[P5C] 


Atkinson 


Radical Chemistry: 










0H*+O3=H0 2 *+0 2 


6.7*10" 14 


cm^ molecule* 1 s -1 


[Rl] 


DeMore 


H0 2 *+0 3 =0H*+20 2 


2 .o*io - 15 


cm3 molecule" 1 s'* 


[R2] 


DeMore 


HC0*-+O2=H02*4C0 


5.6*10" 12 


cm3 molecule" 1 s" 1 


[R9] 


Anatasi 


oh*+h 2 o 2 =h 2 o+ho 2 * 


1.7*10" 12 


cm3 molecule -1 s -1 


[R13] 


Seinfeld 


H*+0 2 =H02* 


1. 2 *10" 12 


cm3 molecule" 1 s" 1 


[R18] 


Anatasi 


ho 2 *+ho 2 *=h 2 o 2 +o 2 


1.6*10" 12 


cm3 molecule" 1 s* 1 c 


[R20] 


Seinfeld 


N0+H02*=N02+0H* 


8.3*10' 12 


cm3 molecule" 1 s" 1 


[R22] 


Seinfeld 


no 2 +ho 2 *=hono+o 2 


1.4*10" 14 


cm3 molecule" 1 s" 1 


[R23] 


Stockwell 


co+ho 2 *=co 2 +oh* 


1.9*10- 3 2 


cm3 molecule" 1 s" 1 


[R24] 


Atkinson 



c Assumed [HO2] to be low, considered R20 to be negligible. 

Use Steady State Assumption: 

[H02*]ss=(P5C+R1+R9+R13+R18)/(KR2*03+KR22*N0+KR23*N02-KR24*C0) 



Carbon Monoxide Reactions (CO): 
Photolysis: 



HCHO+hv -*• CO+H? 


>j>300nm 


3.86E-5 s" 1 


[PI B] 


Atkinson 


Radical Chemistry: 


HCO*+O 2 =HQ 2 *4C0 


5.6*10 -12 


cm 3 


molecule" 1 s" 1 


[R9] 


Anatasi 


HCHO+OH*= H*+C0+H 2 0 


1.1*10" 11 


cm 3 


molecule" 1 s" 1 


[RIO] 


Seinfeld 


COOH*=C07+H* 


2. 2 *10 -13 


cm 3 


molecule -1 s" 1 


[R13] 


DeMore 


COt-HCh*= C02+0H* 


1.9*10" 32 


cm 3 


molecule" 1 s" 1 


[R24] 


Atkinson 



d(CO]/dt=PlB+R9+R10-R15-R24 



Atomic Hydrogen Reactions (H*): 
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Photolysis: 



HCHO+hv -*H*+HCO* 


X<300nm 


1.30E-5S - 1 


[PI A] 


Atkinson 


H 2 O 2 +I 1 V -h*+ho 2 * 


k<324nm 


0, (4> = 0) 


ffsq 


Atkinson 


H202+hv -2H*+02 


X<213nm 


0, (4> = 0) 


[PSD] 


Atkinson 


H 2 0+hv — H*+OH* 


k<239nm 


1.09E-5 s' 1 


[P7B] 


Atkinson 


HONO+hv — H*+N02 


X<367nm 


0 , (<t> = 0) 


[P50B] 


Atkinson 


Radical Chemistry: 


HCHOOH*=H*4COH20 


l.mO* 11 cm 3 


molecule' 1 s' 1 


[R10] 


Seinfeld 


CO+OH*=C02+H* 


2.9*10* 13 cm 3 


molecule* 1 s' 1 


[R1S] 


DeMore 


N+OH*=NO+H* 


4.9*10- u cm 3 


molecule' 1 s' 1 


[R10] 


Atkinson 


H*+0 2 = HO 2 * 


1.2*10* 12 cm 3 


molecule" 1 s' 1 


[R18] 


Anatasi 


H*+03 = OH*+ 0 2 


2.9*10' 1 1 cm 3 


molecule* 1 s* 1 


[R19] 


Anatasi 



d[H*]/dt=PlA+P5C+P5B+P7B+P50B+R10+R15+R16-R18-R19 



Formyl Radical Reactions (HCO*): 
Photolysis: 



HCHO+hv — H*+HC0* 


X<300nm 1.30E-5s * 1 


[PI A] 


Atkinson 


Radical Chemistry: 


HCO*+0 2 = H02*4<X) 


5.6*10* 12 cm 3 molecu v 1 s* 1 


[R9] 


Anatasi 


HCHO+OH*= H 2 O+HCO* 


1.0*l0* n cm 3 molecule* 1 s* 1 


[R17] 


Atkinson 



d[HCO*]/dt=PlA-R9+R17 
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Atmospheric Gases: [Seinfeld, p. 8] 
Gas Avg. Concentration rppml 



Ar 

Ne 

Kr 

Xe 

N 2 

Oz 

CH4 

CO 2 

CD 

H? 

N 2 0 

S0 2 

nh 3 

NO+NO 2 

O 3 

HNO 3 

CH 

HO 2 

Os 

H?0? 

H?0 

He 


9340 

18 

1.1 

0.09 

780,840 

209,460 

1.65 

332 

0.05-0.2 

0.58 

0.33 

10-5-10-4 

10-4-10-3 

10-6-10-2 

lO-2-io-i 

10-5-10-3 
0.406*10‘ 5 ppb 
0.406*10' 2 ppb 
30 ppb 
10 ppb 
22,550 
5.2 
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